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LONG-TERM  GOALS 

We  have  the  long  term  goals  of  understanding  several  aspects  of  passive  and  active  sonar  performance 
in  the  Philippine  Sea  when  there  is  a  robust  environmental  characterization  of  the  sea  floor  and  water 
column  plus  accurate  source/receiver  positions.  This  can  lead  to  better  signal  processing  where  the 
structure  and  coherence  can  be  exploited. 

OBJECTIVES 

The  performance  of  the  beamforming  of  a  sonar  can  be  separated  into  the  gain  against  noise,  or  NG  and 
the  signal  gain  degradation,  or  SGD.  These  two  quantities  are  well  understood  for  horizontal  line  arrays 
or  arrays  without  much  vertical  extent;  however,  the  same  cannot  be  said  when  using  large  vertical 
arrays.  The  vertical  noise  structure  remains  and  ongoing  effort  for  VLA’s  Here  we  are  concerned  with 
sources  of  the  SGD.  Past  efforts  in  ocean  acoustics  concerning  the  vertical  structure  were  termed 
matched  field  processing  (MFP)  or  coherent  multipath  recombination  (CMR).  (1)  Both  of  these 
approached  have  not  fulfilled  their  potential  for  a  number  of  reasons  including  unrealistic  assumptions 
about  the  coherence  of  the  vertical  field  and  especially  the  multipath  and/or  mode  coherence. 

The  Philippine  Sea  experiments  deployed  two  large  vertical  arrays,  or  DVLA’s  (Deep  Vertical  Line 
Arrays).  The  first  experiment  in  2009  had  two  sections  at  the  upper  and  lower  turning  points  while  the 
second  in  2010  had  140  sensors  which  spanned  the  5600  meter  water  column.  (4)  The  objective  here  is 
to  determine  use  these  data  to  address  these  coherence  issues.  The  first  issue  is  assess  the  impact  of  an 
internal  wave  field  on  the  randomizing  of  the  multipaths.  This  issue  is  controversial.  A  signal  from  a 
point  source  starts  as  a  fully  coherent  signal  field  before  any  random  medial  effects  can  accumulate. 

The  issue  then  is  at  what  range  does  the  media  lead  to  uncorrelated  multipath  which  is  the  favorite 
model  for  many  in  acoustical  oceanography.  Nevertheless,  an  uncorrelated  vertical  field  essentially 
dooms  MFP  or  CMR  processing.  If  the  field  is  coherent  for  a  few  convergence  zones,  this  is 
consequential  to  potential  of  most  low  frequency  sonar  systems  with  vertical  apertures. 
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APPROACH 


The  DVLA  sections  for  the  first  experiment  were  close  enough  for  beamforming  in  the  vertical; 
however,  those  for  the  second  had  40  meter  separation  which  leads  to  an  aliasing  angle  of 
approximately  8  degrees.  Because  of  all  the  calibrations  needed  to  interpret  the  data  correctly  we  have 
use  simulations  methods  to  guide  our  observations.  This  has  been  done  by  using  an  internal  wave  code 
implementing  the  Garrett-  Munk  model  for  randomizing  the  media  and  then  propagating  the  RAM  wide 
angle  parabolic  equation  code  through  these  media.  1 

The  approach  is  to  simulate  the  signals  observed  on  a  VLA  generated  by  an  axial  source.  The  VLA 
receiver  model  has  a  1/2  wavelength  spacing  at  250  Hz,  which  is  the  center  frequency  for  many  of  the 
Philippine  Sea  transmissions  and  extends  throughout  the  watercolumn  to  the  ocean  bottom.  2  We  then 
consider  a  sequence  of  subapertures  with  sixty  four  sensors  which  leads  to  each  one  have  an 
approximately  two  degree  resolution.  Figure  1  illustrates  the  array  geometry  and  a  sample  PE  solution. 
The  first  step  is  first  to  identify  the  path  structure  at  each  range.  For  this  apply  both  conventional  and 
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Figure  1:  Left:  Concept  of  overlapped  apertures  on  a  VLA.  Right:  Example  of  RAM  solution  w/ 

Garrett-Munk  perturbation 

adaptive  array  algorithms  to  each  subaperture,  or  a  directional  wavenumber  spectrum  spectrum 
parameterized  by  the  depth  of  each  subaperture.  For  this  we  consider  each  subaperture  to  have  64 
elements,  or  an  approximate  resolution  of  2  degrees.  The  apertures  are  overlapped  by  32  elements,  so 
each  step  in  depth  corresponds  to  an  increment  of  48  meters. 


RESULTS 

Figure  2  illustrates  two  of  these  depth  dependent  spectra  for  ranges  of  120  km,  or  two  convergence 
zones  (CZ’s),  and  150  km,  or  two  and  one  half  CZ’s.  We  label  such  figures  as  a  depth-direction  sine 

1  Both  the  internal  wave  and  RAM  MATLAB  codes  were  provided  by  Matt  Dzieciuh  of  Scripps  Institute  of  Oceanography. 

(3) 

“A  spacing  of  a  1/4  wavelength  would  suffice  because  high  angle  wave  components  were  suppressed  by  an  absorbing 
bottom  boundary. 
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spectrum.  (While  CZ  propagation  is  usually  discussed  in  the  context  of  shallow  sources  and  receivers, 
the  RAM  solution  reflects  a  CZ  like  behavior.)  The  figures  have  encapsulate  the  ray  arrivals  very 
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Figure  2:  Left:  Directional  wave  spectra  vs  subaperture  depth  at  120  km;  Right:  Directional  wave 
spectra  vs  subaperture  depth  at  150  km.  Note  the  subaperture  depth  is  on  the  vertical  and  direction 
sine  on  the  horizontal  and  is  measured  normal  to  the  VLA,  or  the  usual  grazing  angle. 

concisely  resembling  what  is  often  termed  and  angle-depth  diagram.  At  the  120  km  range  we  can 
identify  two  paths  at  the  shallow  subapertures  corresponding  to  signals  refracted  near  the  sea  surface. 

As  we  progress  deeper  there  are  four  paths  corresponding  to  two  additional  signals  for  the  signals 
ducted  along  the  sound  channel  axis  at  1000  m  where  they  were  launched  at  the  source  depth.  Going 
deeper  near  1500  m  we  can  identify  four  paths  up  to  2000  km.  Finally  we  loose  the  high  angle  paths 
and  observe  just  those  signals  turning  at  low  grazing  angles.  The  depth  dependent  spectra  for  150  km, 
or  between  a  CZ  range,  has  a  quite  different  structure.  First,  there  are  no  shallow  rays  near  the  surface. 
Significant  energy  does  not  appear  until  approximately  250  m.  After  this  there  are  just  two  paths  until 
1500  m.  One  is  a  high  angle  diving  path  with  a  negative  D/E,  while  the  other  is  a  packet  of  positive  D/E 
angles  spread  from  grazing  to  nearly  a  direction  sine  of  .2,  or  approximately  10  degrees.  After  this  there 
is  but  one  signal  with  a  negative  D/E  until  the  critical  depth  of  the  profile  when  two  paths  reappear.  A 
very  interesting  analysis  can  be  done  while  observing  the  stepping  out  in  range  as  the  structure  of  the 
spectra  changes.  Space  does  not  allow  including  such  a  figure.  Eventually,  we  expect  that  the 
propagation  will  ’’homogenize”  with  a  mostly  symmetric  depth-direction  sine  spectrum.  Eventually, 
bottom  interaction  ceases  to  absorb  energy.  Then,  the  volume  scattering  introduced  by  the  internal  wave 
perturbations  then  will  lead  to  an  ’’equipartition”  of  all  the  propagating  rays  and/or  modes.  At  these 
ranges  there  is  no  potential  to  exploit  the  vertical  structure.  One  way  of  quantifying  the  coherence  of  the 
signals  both  at  each  subaperture  and  the  among  them  involves  using  the  covariance  matrices  for  the 
array  data  at  eace  subaperture.  Two  measures  using  thess  are  the  mutual  information  and  the  cross 
coherence  among  the  subapertures.  There  is  in  fact  a  one  to  one  relationship  between  the  two  measures. 
The  mutual  information  is  a  quantitative  measure  of  how  much  information  one  signal  vector  provides 
about  another.  It  is  given  by 

/(X!,X2)  =  -E[PxuX2(XhX2)\. 
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For  complex  Gaussian  random  vectors  this  is  given  by 


/(Xi,X2)  =  log2 
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,  where  £  denotes  the  auto-covariances  and  cross-covariances  according  to  the  subscripts  and  the  |  •  |  is 
the  determinant  operator.  By  using  the  properties  of  partitions  of  determinants  and  their  combination 
leads  to  an  expression  of  the  mutual  information  in  terms  of  the  magnitude  of  the  coherence  metric,  or 

/(X1?X2)  =  -/o,?2[l-|p(X1,X2)|2] 


with 

|P(Xi,X2)|2 

which  is  the  multivariate  definition  of  coherence. 
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We  have  applied  this  to  covariances  computed  by  using  the  ensemble  of  the  PE  realizations  on  the 
subapertures.  3  Another  measure  of  coherence  is  the  number  of  degrees  of  freedom  (dofs) for  an 
ensemble  of  signal  vectors  observed  on  an  a  subaperture  as  well  as  for  a  pair  of  subapertures.  For  a 
single  subaperture  a  measure  of  this  is  given  in  terms  of  the  eigenvalues  of  a  covariance  as  well, (5) 
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where  the  A,  are  the  eigenvalues  of  the  matrix  E. 

We  are  applying  all  these  measure  to  assess  the  wavefield  coherence.  In  the  near  future  we  plan  to 
implement  the  two  measures  described.  After  we  understand  the  issues  using  the  simulated  data, 
especially  the  statistical  and  calibration  ones,  we  plan  to  applying  this  to  the  DVLA  data  of  the  Phil  Sea. 


IMPACT/APPLICATIONS 


There  are  a  number  of  ongoing  plans  for  the  use  of  large  vertical  apertures  for  Navy  systems.  Some  of 
these  are  at  the  low  frequencies  used  in  the  Phil  Sea  experiments,  while  other  are  at  mid  frequencies  of  a 
few  kHz.  The  aperture  scales  mostly  scales  as  the  inverse  of  the  frequency;  nevertheless,  the  coherence 
of  the  multipaths  is  an  important  factor  in  the  design  of  signal  processing  for  VLA’s.  All  current 
algorithms  assume  incoherent  paths,  similar  to  the  distribution  in  the  horizontal  originating  from 
uncorrelated  sources.  An  important  implication  is  that  there  may  be  less  SGD  and  higher  NG  once  we 
understand  the  statistics  of  the  multipath  coherences. 


3  There  is  a  significant  statistical  issue  for  computing  these  covariances  beyond  the  discussion  here.  It  is  often  suggested 
that  one  needs  sample  support  on  the  order  of  3X  the  number  of  sensors  which  implies  approximately  200  realizations  for 
statistical  stability  of  the  low  order  eigenvalues.  This  is  very  much  an  issue  for  a  determinant  since  it  is  given  by  the  product 
of  the  eigenvalues. 
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